Rabi oscillations of 2D electrons under ultrafast intersubband excitation 



o 
o 

(N 



^3 

o 

CD 



CO 



I 

C 

O 

o 



> 

(N 

m 

in 
o 
m 
O 



13 

o 
o 



X 
S3 



D. McPeakeQ F. T. Vasko, and E. P. O'Reilly 
NMRC, University College, Lee Mailings, Prospect Row, Cork, Ireland 
(Dated: February 2, 2008) 

We investigate coherent nonlinear dynamics of 2D electrons under ultrafast intersubband excita- 
tion by mid-IR pulses. We include the effects of relaxation and dephasing, both homogeneous and 
inhomogeneous, as well as detuning within a non-Markovian equation to obtain temporal population 
redistributions. We show how, using a cross-correlation method, the effects of Rabi oscillations may 
be detected in this system, and briefly discuss other detection methods. 
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The coherent dynamics and manipulation of two-level 
systems under ultrafast excitation dates back over three 
decades, where the majority of focus has been on atomic 
ensembles 0. During recent years, much attention 
has been given to interband coherent effects in low- 
dimensional heterostructures 0- In particular, investi- 
gations of single quantum dots under interband ultrafast 
excitation |3|, where dephasing times can be hundreds 
of picoseconds, has seen the experimental observation of 
Rabi oscillations, providing the basis for possible readout 
devices for excitonic quantum gates Q . Previously, inter- 
band Rabi oscillations under near-IR excitation had been 
observed in quantum well (QW) structures, albeit more 
highly suppressed by shorter dephasing times 0, . Fol- 
lowing these interband investigations we investigate the 
possibility of Rabi oscillations between n = 1 and n = 2 
conduction subbands in a QW under ultrafast mid-IR ex- 
citation. Advances in mid-IR femtosecond spectroscopic 
techniques 7] have culminated with recent experimen- 
tal results y| suggesting that 'partial Rabi flops' may be 
possible under relatively low pulse energies. Theoreti- 
cal investigations of intersubband systems are therefore 
timely. 

The occurrence of subband Rabi oscillations under rel- 
atively low pump energies in low-doped QWs is made 
possible due both to the suppression of the LO-phonon 
emission rate (assuming the interlevel energy is much 
greater than the optical phonon energy) and also the 
reduced efficiency of electron-electron scattering at low 
densities (~ 10 9 — 10 10 cm -2 ). The inclusion of the 
role played by dephasing mechanisms is crucial for any 
description of coherent intersubband optics. Whilst a 
complete treatment requires the full rigour of many- 
body calculations, we show in this paper that treating 
the different damping terms phenomonologically within 
a non-Markovian quantum kinetic scheme gives clear 
indications of the interplay between the various time- 
scales involved, and provides insight into the possiblity 
of coherent manipulation. Experimental observations 
9] show that carrier-carrier dephasing rates are effec- 
tively independent of carrier density at low concentra- 



tion (~ 10 10 cm -2 ). We therefore restrict our study to 
the low-density regime where the assumption of constant 
dephasing is approximately valid. 

After introducing the framework of the two subband 
model, we consider the balance of population redistribu- 
tions to derive an integro-differential equation governing 
population dynamics under the influence of an ultrafast 
pulse. By considering the ideal limit, where dephasing 
and relaxation are neglected, we show that our model 
reduces to analytic expressions which essentially recover 
the area theorem ^(j . The effects of the various damping 
mechanisms are reintroduced in line with recent experi- 
mental measurements to investigate the degree of coher- 
ent behaviour which can realistically be observed. We ex- 
tend our calculations to analyse a detection scheme based 
on cross-correlation absorption under two-pulse excita- 
tion, and briefly discuss a scheme based on THz emission 
from non-symmetric QWs. Such methods have previ- 
ously been implemented for the near-IR interband pump 
case 0, . We conclude with a brief discussion on the 
approximations we have made. 

We describe the coherent dynamics of the system us- 
ing a density matrix approach, evaluating the balance 
equation for an interlevel redistribution of population 
under a resonant intersubband pump. We consider the 
response of the system due to the interaction with an 
electric field E t = [i?exp(— iurt) + c.c.jwt perpendicular 
to the QW plane. The dipole perturbation operator is: 
5ht exp (— iut) + H.c, with Sht = (ieE/uj)v±Wt, where 
v± is the intersubband velocity operator and Wt is the 
form-factor of the pulse. We write the high-frequency 
component of the density matrix, Spt exp {—iut) + H.c, 
in the form |l2j : 
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where t%, the dephasing time, describes the decay of the 
induced coherence, and h is the period-averaged Hamil- 
tonian of the QW. The density matrix averaged over the 
period 2-7T /lj, p t , obeys the quantum kinetic equation 
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where I sc is the collision integral. The time-dependent 
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generation rate, Gt, is given by: 
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with non-Markovian dependence on p t . In the resonant 
approximation, we consider two states \jp), j = 1,2, 
which we take to have parabolic dispersion laws, Sj + 
p 2 /2m. In this case one can neglect the off-diagonal terms 
in the density matrix under the condition: h/e 2 \T <C 1, 
where £21 is the interlevel energy and r is the character- 
istic relaxation time. According to particle conservation, 
we have that: G\ v t — —G 2p t = G p t, and the generation 
rate can be written in the form: 

\ Tkf ) m J dt ' w t+v<t>t> (/2 P t+f - Apt+t') • 

(4) 

Here fj pt = (jp\pt\jp) is the population of the jth state, 
(2|0j_|l) = (1|ujJ2)* = iv± is the interband velocity, 
Alu = Lu~e 2 \/h the detuning frequency, and the kernel is 
given by <pt = cos(Aui) exp (t/Vjj). Taking into account 
the inhomogeneous broadening of intersubband transi- 
tions we consider |l3j an in-plane nonuniform intersub- 
band energy e 2 i ( x ) = £21 in the x-dependent gener- 
ation rate G pX f By averaging over the 2D plane, the de- 
tuning term in the generation rate becomes: (cos Aoj x t). 
We therefore replace the kernel <fi t in Eq. (4) with 
$ t = <j) t exp[— (-yt/^/2H) 2 ] where 7 = yj (Se x Se x ) is the 
energy broadening due to nonuniformity. 

Within this framework, we can now define the concen- 
tration in the jth level as nj t = (2/ L 2 )J2 P fjp.u L 2 is 
the normalisation area and the factor 2 is due to spin. 
The balance equations have the form: 
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where G t = (2/L 2 )^2 p {jp\G t \jp) and (dn 3 /dt) sc = 
(2/ 1 L 2 )^ p {jp\I sc \jp) describes interlevel relaxation. 
Due to particle conservation, n\t + n 2t = "B2D where 
ns2D is the total 2D concentration, so that (dni/di) + 
(dri2/dt) 8C = 0. Defining the population redistribution 
An t — nit — n^t, and describing the interlevel relaxation 
by the relaxation time n, we transform Eq. (5) into 



dAn t 
dt 
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n 



(6) 



where the generation rate Gt is expressed through An t 
according to Eq. (4). We solve Eq. JHJ subject to the 
initial condition Ant^-00 = n>&2D- Upon substitution of 
Eq. (4) into (6) we obtain the non-Markovian integro- 
differential equation 



dAnt 
dt 
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where we have introduced the phototransition frequency 
v r = (2ev±E/hw) 2 T p with r p the pulse duration. Factors 
which suppress the coherent response are described by 
the parameters (ti^) -1 and j/h, as well as detuning of 
the pulse from resonance Alu. 

Before we resort to numerical solution of Eq. (7) , we 
consider analytic solutions in the limiting case Tx,i —* 00, 
7^0, and Alu — > 0, ( i.e. $ t — > 1). Eq. (7) can then be 
transformed into the second order differential equation 



d 2 An t 1 dwt dAn t 
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wt dt dt 
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to which we apply an additional boundary condition: 
dAnt/dij t ^_ QO = 0. The population dynamics in 
this regime can be described by the expression 



t/T 



Ant = cos yJv r T v 



dzw. 



(9) 



which we illustrate in Fig. 1. As expected, population 
redistributions are solely determined by the area of the 
incident pulse: A p = ^Jv r r p j_ oQ dzw z . In Fig. la we 
plot the solution as a function of time for both secant 
and Gaussian pulses each with duration t p (FWHM) . As 
illustrated in panels (i) and (ii), which correspond to se- 
cant pulses with A p = tt and 2.5tt respectively, popula- 
tion can be distributed among the two levels according 
to the area theorem. For fixed values of v r and r p , the 
area under the Gaussian pulse is slightly less than that 
under the secant pulse. The final redistribution of popu- 
lation, shown in Fig. lb as a function of the secant pulse 
area, Tr^/v r T p / 'a s , where a s = 2 sech _1 [0.5], confirms this 
point. 

The condition for population inversion, A p — it can 
be obtained by a secant pulse, with duration t p =100 
fs, having a free space intensity per unit area of 5.22 
MW/cm 2 . In order to achieve maximum coupling of the 
transverse field to the intersubband dipole, we consider 
a 45° prism integrated onto the sample. With this ge- 
ometry, approximately 33% of the incident power can be 
coupled to the transition. We assume parameters for a 
typical GaAs QW, of width 85 A, corresponding to an in- 
terlevel energy of Too ~100 meV. Both sets of parameters 
are typical of recent experiments y\. 

To study the effects of the various damping processes, 
we solve Eq. (7) numerically using a form of Picard it- 
eration. We take as example a secant 1.67r-pulse, and 
consider the effects of the parameters, t±, t 2i Aw -1 , and 
v27i/7, setting each in turn to t p with the other three set 
to the limiting values as outlined above. The results are 
shown in Fig. 2. LO-phonon emission (n) causes slow re- 
laxation of the population back to the ground state after 
the pulse has passed. In contrast, the effect of a large de- 
tuning as well as broadenings, both inhomogeneous and 
homogeneous (y^h/j, r 2 ), is to cause suppression of the 
Rabi flop during the pulse. Clearly the potential for co- 
herent manipulation depends on the interplay between 
these processes. 
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FIG. 1: (a) Temporal population redistribution, Ant, for the 
ideal case (ti — > oo, $t — > 1) due to secant (thick solid curve) 
and Gaussian (dashed curve) pulses of duration t p (FWHM). 
The pulse profiles are shown as the thin solid (secant) and 
dotted curves (Gaussian). Panels (i) and (ii) correspond to 
pulse areas (for secant pulse) A p = it and 2.5n respectively, 
(b) The area theorem is recovered as final distribution is plot- 
ted against secant pulse area (1.1 93^/^ r r p ) for secant (solid 
curve) and Gaussian (dashed curve) profiles. 



To include realistic estimates of these parameters, we 
take experimentally determined values of coherence times 
for a typical sample [8J, |9| . Four- wave- mixing measure- 
ments provide a value for the homogeneous broadening 
corresponding to an electron-electron scattering time of 
the order 320 fs. LO-phonon emission times are esti- 
mated to be between 1-4 ps. As a conservative estimate 
we take n = 1 ps and calculate the system response to 
a 100 fs (FWHM) pulse tuned to the 1 -> 2 transition at 
2?= 100 meV, and allow for a detuning of up to 12 meV. 
In Figs. 3a, b we show the population redistributions due 
to these parameters. Although complete inversions are 
suppressed, we estimate that up to 70-80% of the initial 
distribution can be coherently excited to the upper level 
with a resonant 7r-pulse: Fig. 3a (i). Effects of pulse 
area can still be seen in the population 1 ps after the 
maximum of the pulse intensity. However the question 
remains as to whether these density flops can be detected 
experimentally. We now consider two possible scenarios. 

Following the interband example p|, we calculate 
cross-correlation functions comparing changes in pulse 
modification due to transmission through unexcited 
and weakly excited samples. With the induced cur- 
rent density written in the standard form: J t — 
(2e/i 2 ) £,j,p v 3 y (j'p\Sp t \jp), we can use Eq. (1) to ex- 
press the induced current through An t . The absorbed 
power, Vt — JtEt where the over line means the average 




FIG. 2: The suppression of the ideal Rabi flop (solid curve) 
under 1.67r-pulse excitation (duration r v (FWHM)) due 
to: detuning (dot-dashed curve), inhomogeneous broaden- 
ing (dashed curve), homogeneous broadening (dotted curve) 
and LO-phonon emission (dot-double-dashed curve). Each 
parameter is varied in turn with the rest set to ideal values. 
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FIG. 3: (a) Temporal population redistributions between sub- 
bands levels 1 and 2 due to excitation by a 100 fs (FWHM) 
secant pulse of area (i) n, (ii) 3tv. Decay parameters are taken 
to be T2 = 320 fs, ri = 1 ps, corresponding to a two-level sys- 
tem homogeneously broadened by 4 meV. Also shown is the 
effect of detunings of 6 meV (dot-dashed curves) and 12 meV 
(dashed curves). The fraction of population excited at 1 ps 
after the maximum pulse intensity is shown as a function of 
A p in (b). 



over the period, is written in the form: 

Vt = - — r^-^wt / dt'wfe 1 ,T2 cosAujt'An t+t <. 
riuJ J-oo 

(10) 

The transmitted intensity of a single pulse is thus given 
by: I s = I t — Vt, where It oc \E t \ 2 . Id, the transmitted 
intensity of an identical pulse after weak excitation of the 
sample by a prepulse, is similarily extracted from Eqs. 
(7)(10) using the two-pulse form-factor, w t — > wf re + 
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FIG. 4: Differential cross-correlations, shifted vertically for 
clarity, for different pulse areas (numbers to the left of peaks) 
as a function of reference time delay, Tfi re /. Inset: Temporal 
absorption (dotted curve) and density variation (solid curve) 
for the pulse with area A p = 1.25n. 

(3wt where /3 = E/E pre is the ratio of the maximum of 
the electric field of the pulse and prepulse. The cross- 
correlation is thus: 

hxc(TKref) OC J dt [I d (t) ~ I s (t)} I r (t ~ T Qre f) , (11) 

where I r (t—TQ re f) is the intensity of a reference pulse (we 
take I r (t — TR re f) oc \E t \ 2 ). The interval between pulse 
and prepulse is kept fixed at 600 fs allowing for no coher- 
ent interaction between pulses. The prepulse area is fixed 
at A^ re ^ = 0.257T. The intensity (area) of the pulse and 
hence reference pulse is varied and the cross-correlation 
(11) calculated. The results are shown in Fig. 4. The ap- 
pearance of zero crossings in the cross-correlation, which 
move to shorter time-delays with increasing pulse inten- 
sity, is exactly the behaviour expected due to Rabi flops. 
These crossings which occur due to the dependence of the 
cross-correlation signal on the temporal derivative of the 
excitation density during the pulse, show that the system 
has been driven through a density maximum - a partial 
Rabi flop. 

In addition to this cross-correlation scheme, we envis- 
age detection based on the charge oscillations of a Rabi 



flopping population in a non-symmetric QW with dipole 
moments d\ and di- The temporal dependence of the 
induced dipole moment, ni t di+ri2td2 oc (1 — dz/dijArit, 
connects the temporal population dynamics, An t , with 
the induced oscillations of the electric field. The response 
can be measured by the method employed in |£|, as well 
as by the emission of THz radiation. We note that the 
experimental method [ll| used for the detection of THz 
radiation under interband excitation was carried out with 
carrier densities comparable to those we have considered. 

Our calculations are based on a few assumptions. Since 
the Coloumb renormalisation of the intersubband ab- 
sorption (due to depolarisation and exchange effects) 
is weak in low-doped QWs |9, [l3j> we have adopted 
a single-particle approach. All homogeneous relaxation 
and dephasing mechanisms including electron-electron 
scattering and LO-phonon emission have been considered 
phenomenologically with inhomogeneous broadening in- 
cluded in the long-scale limit. We expect that micro- 
scopic considerations based on many-body calculations 
of the various damping processes, taking into account in- 
plane inhomogeneities, may slightly modify the nonlinear 
response. However, we do not expect large departures 
from this simple model for the low-density regime we 
have considered. Due to the prism geometry we have as- 
summed, the finite length of the pulse, L, must be taken 
into account. For typical experimental parameters we 
estimate that L w 15 um, thus setting a limit on the 
sample area, and hence amplitude of response we can de- 
scribe. However, we note that the temporal dynamics 
we have described are not changed as L typically exceeds 
the characteristic kinetic lengths, vt\^ and v/v r , where 
v is the characteristic in-plane velocity. 

In conclusion, we have developed a simple description 
of Rabi oscillations between subband levels in low-doped 
QWs. We have considered the population dynamics of 
a two-level system under ultrafast mid-IR radiation in 
a non-Markovian quantum kinetic model. Our calcula- 
tions, based on realistic pulse and material parameters, 
show that coherent manipulation of the order of 10 6 elec- 
trons (depending on area of excitation) should be possible 
under relatively low-power excitations. 
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